Optical polarization and near IR photometry of the 
proto-planetary nebula Hen 3-1475^ 
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ABSTRACT 

We present BVRI CCD aperture polarization and near-infrared photometry of the proto- 
planetary nebula Hen 3-1475. Its intrinsic polarization is high and shows a strong spectral 
dependence. The position angles in all bands are perpendicular to the axis of the observed 
bipolar structure. A Monte Carlo code is used to model the intrinsic polarization of Hen 3- 
1475. Using disk dimensions and other constraints suggested by previous works, we are able to 
reproduce the observations with an optically thick disk composed by grains with a power-law size 
distribution ranging from 0.06 to 0.22 /zm. We also reliably estimate the foreground polarization 
from hundreds of stars contained in the CCD images. It is parallel to the intrinsic polarization of 
Hen 3-1475. Possible implications of this result are discussed. From IR observations, we estimate 
a interstellar reddening, Ay, of about 3.2. 

Subject headings: stars: individual (Hen 3-1475) — planetary nebulae: individual (Hen 3-1475) — 
circumstellar matter — polarization — radiative transfer 



1. Introduction 

Hen 3-1475 (IRAS 17423-1755) is identified as 
a proto-planetary nebula (PPN), i.e., an object 
between the end of the asymptotic giant branch 
(AGB) and the planetary nebula (PN) phases. 
The classification and physical properties of Hen 3- 
1475 are supported by a broad base of observa- 
tional data: radio emission (Knapp et al. 1995); 
OH masers (te Lintel Hekkert 1991; Zijlstra et al. 
2001); optical imaging and spectroscopy (Riera 
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et al. 1995; Bobrowski et al. 1995; Borkowski, 
Blondin, & Harrington 1997; Borkowski & Har- 
rington 2001). Probably the strongest evidence in 
favor of its evolutionary stage is the high N abun- 
dance, typical of evolved objects, specifically of 
Type I PN (Riera et al. 1995). The radio contin- 
uum indicates the presence of a compact ionized 
region surrounding the central B star with an es- 
timated dynamical age of only 15 years (Knapp 
et al. 1995. Recently, Borkowski & Harrington 
(2001) have obtained images and detailed spec- 
troscopic observations and, among other results, 
they have derived a distance of 8 kpc, which im- 
plies a very high luminosity (25 000 L©) for a AGB 
star at the Galactic Bulge. Analyzing Ha spec- 
tra, Sanchez Contreras & Sahai (2001) could iden- 
tify two winds, one of which is probably associ- 
ated with the post-AGB wind before its interac- 
tion with the previous AGB wind. 

High-resolution optical images of Hen 3-1475 
show a very complex bipolar morphology (Borkowski 
et al. 1997; Bobrowski et al. 1995; Riera et al. 
1995). The central region is dense, probably a 



1 



dust disk, which may be responsible for the OH 
maser emission as well as for the Call and Fell 
lines observed in its near-infrared spectrum. This 
structure has an angular size of about 2". The 
high resolution image obtained by Borkowski et 
al. (1997) shows also dark patches probably asso- 
ciated with denser clumps of material in the disk. 
A complex structure of highly coUimated optical 
jets with pairs of point-symmetrical shock-excited 
knots produces a characteristic S-shapcd sequence, 
which is about 15" in size oriented along PA 135°. 

In the rapid evolutionary phase of PPN, dra- 
matic changes in geometrical and physical prop- 
erties of its central star and circumstellar mate- 
rial take place. An example is the development 
of non-symmetric geometries (eg., bipolar nebu- 
lae) in post-AGB objects. The interaction be- 
tween an equatorial enhanced AGE wind and that 
from an early PN phase was suggested as a mech- 
anism to produce the variety of geometries seen 
in PPNe and PNe (Kahn & West 1985). The sta- 
tus of this and other models is reviewed by Frank 
(2000). This enhanced equatorial wind - which 
we will call disk - can also play a role in colli- 
mating the jets seen in many objects. Recently, 
Frank, Lery, & Blackman (2002) have shown that 
a magnetic-centrifugal launching model can ex- 
plain winds with velocities as high as those ob- 
served in Hen 3-1475. Therefore, a good repre- 
sentation of the disk properties is crucial to un- 
derstand the changes occurring during the PPN 
phase. 

Polarimetry is an useful tool in the study of 
PNe and PPNe because it is sensitive to the cir- 
cumstellar geometry and composition. The polar- 
ization measured in PPNe is usually high due to 
the scattering in asymmetric circumstellar envi- 
ronments. .Johnson & Jones (1991) have shown 
that, from red giants to PNe, the evolutionary 
stage of PPN is the one in which the polarization 
tends to be highest and that the transition from a 
spherically symmetric envelope (typical of red gi- 
ants) to a non-spherical envelope seems to occur 
in early stages of the AGB phase. The wavelength 
dependence of the polarization in PPNe has vari- 
ous shapes and some objects present polarization- 
angle rotation (Trammell, Dinerstein, & Goodrich 
1994). Many extended objects were observed us- 
ing imaging polarimetry (e.g., Scarrot & Scarrot 
1995; Gledhill et al. 2001). Far from the central 



source, the light is highly polarized and exhibits 
the centro-symmetric pattern typical of reflection 
nebulae. However, the polarization near the cen- 
tral object tends to be parallel to the optically 
thick disk originated in the AGB phase. This 
pattern is also observed in young stellar objects 
with bipolar geometries (see references in Bastien 
& Menard 1988 - Table 1). 

In spite of the importance of grain size in 
the study of the temperature profile and, conse- 
quently, of the infrared (IR) luminosity (Stasinska 
& Szczerba 1999, 2001), few works have focused 
on the determination of the grain size distribution 
in PPNe and PNe. An example of model that con- 
siders the grain sizes as free parameters is found 
in Meixner et al. (2002). That work, however, 
modeled the observed spectral energy distribution 
(SED) alone. Recent results (Carciofi, Bjorkman, 
& Magalhaes 2002; Li & Lunine 2002) show that 
the SED alone does not constrain the grain sizes. 
The spectral dependence of the polarization helps 
to restrict the models to the dust size distribution 
(Carciofi, Magalhaes, & Bjorkman 2002). Heckert 
& Smith (1988) modeled the polarization of OH 
0739-14 considering single scattering in an enve- 
lope composed by grains of a single size. Johnson 
& Jones (1991) have modeled the observed polar- 
ization of some PPNe, but they used a fixed MRN 
size distribution ranging from 0.005 to 0.10 fxm. 
To our knowledge, there is no previous attempt to 
describe the polarization data of a PPN with a free 
distribution of grain sizes and multiple scattering. 

In this work, we present data on broad band 
CCD aperture polarimetry as well as near-infrared 
(NIR) photometry of Hen 3-1475 (Section 2). In 
Section 3 the NIR photometry is analyzed in order 
to provide information on the interstellar (IS) and 
intrinsic reddening. The polarimetric observations 
are discussed and modeled in Section 4 providing 
some information about the Hen 3-1475 dusty en- 
velope. In the last section, the main conclusions 
are listed. 

2. Observations 

2.1. Optical Polcirization 

The polarimetric observations have been done 
with the 0.60-m Boiler & Chivens telescope at Ob- 
servatorio do Pico dos Dias, Brazil, operated by 
the Laboratorio Nacional de Astrofisica (LNA), 
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Brazil, using a CCD camera modified by the po- 
larimetric module described in Magalhaes et al. 
(1996). The CCD array used was a SITe back- 
illuminatcd, 1024 x 1024 pixels. The above tele- 
scope and instrumentation gives a field-of-view of 
10'. 5 X 10'. 5. The polarization in the V band was 
measured on 1999 April 11 and July 29; the ob- 
servations in the (RI)c bands were also obtained 
on 1999 July 29 and in the B band, on 1999 July 
30. In each night, we observed polarimetric stan- 
dards stars (Serkowski, Mathewson, & Ford 1975; 
Bastien et al. 1988; Turnshek et al. 1990) in order 
to calibrate the system and estimate the instru- 
mental polarization. The measured values of the 
unpolarized standard stars were consistent with 
zero within the errors: consequently no instrumen- 
tal correction was applied. Measurements using a 
Glan filter were also performed to estimate the ef- 
ficiency of the instrument. They indicate that no 
instrumental correction is needed. 

The images have been reduced following the 
standard steps of differential photometry using the 
IRAF^ facility. Counts were used to calculate the 
polarization using the method described in Ma- 
galhaes, Benedetti, & Roland (1984). The mea- 
sured polarization of Hen 3-1475 is shown in Table 
1. The aperture radius that minimizes photome- 
try errors is 4 pixels (corresponding to 2.44" in the 
sky). It was kept constant for all bands. In the 
Discussion and Conclusion sections we will have 
in mind that the polarization measurements refer 
to an extended region in the vicinity of the central 
object in Hen 3-1475. The multicolor polarimetry 
was obtained in nights of poor photometric condi- 
tions with high sky counts. Although this does not 
have important consequences on polarization re- 
sults, which are based on a differential technique, 
this precludes us from detecting the nebula sur- 
rounding the central object. 

2.2. IR Photometry 

Hen 3-1475 was observed with the CamlV NIR 
imager in the 1.6-m telescope of LNA on 1999 
July 29. The detector is a HAWAII array, with 
0.24" /pixel spatial resolution, covering 4'x 4'. J 



'IRAF is distributed by National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc., under contract with 
the National Science Foundation. 



(1.25/im) and H (1.65/im) band images were ob- 
tained. This region is already available in the Two 
Micron All Sky Survey (2MASS) survey, which we 
used to calibrate our photometry. Our combined 
images have an effective exposure time of 210 s in 
J and 40 s in H and provide S/N « 10 for H=17. 
This is deeper than the 2MASS photometry at 
the same S/N, also with a better spatial resolu- 
tion, allowing us to discuss the gross features of 
the reddening in the line-of-sight to Hen 3-1475. 
The J and H magnitudes for the target were ob- 
tained from short exposure time images and are 
J = 9.75 ± 0.01 and H = 8.44 ± 0.01. 

3. Inter- and Circumstellar Reddening of 
Hen 3-1475 

We have obtained the magnitudes of 213 ob- 
jects around Hen 3-1475. An examination of the 
color-magnitude diagram (CMD) H x (J-H) al- 
lows us to put limits on the total reddening to 
this line of sight, since the distribution of unred- 
dened stars is relatively localized. Hen 3-1475 has 
I = 9.8°, b = -1-5.3°, corresponding to a Hne-of- 
sight that intercepts a relatively dense part of the 
bulge of the Galaxy. The objects we observe with 
H < 14 and (J-H) « -1-0.8 probably are red giants 
in the bulge. If they were not subject to reddening, 
the distribution of their color indices should have 
a maximum around J-H « -1-0.5. This can be ob- 
tained by a simple implementation of the model of 
Wainscoat et al. (1992). Using the Cardelli, Clay- 
ton. & Mathis (1989) parametrization of the IS 
extinction curve, we obtain E(J-H) = 0.09 A(V), 
for R = 3.1. So we conclude that the reddening 
to the red giants of the bulge is A(V) w 3.2 mag 
in this line-of-sight. If R = 5, A(V) would change 
to 2.8 mag. These values are a bit larger than the 
previous estimate of Riera et al. (1995), A(V) < 
2.0 ± 0.4 mag. 

The total reddening to Hen 3-1475 presents a 
much more complex problem, since the object is 
shrouded by circumstellar material. We assume 
an spectral type of B3 (Knapp et al. 1995). The 
object position in the CMD is not consistent with 
that of a dwarf, rather, a supergiant star is needed. 
Considering the color of a B3I (Tokimaga 2000), 
we can estimate a circumstellar E(J-H) of 0.84 af- 
ter subtracting the interstellar value. Converting 
the derived E(J-H) to A(V) is not straightforward. 
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A first guess could be made using a standard inter- 
stellar extinction law: this gives an intrinsic A(V) 
= 9.1 (7.9), if R = 3.1 (5.0). However, the spectral 
energy distribution of Hen 3-1475 may be quite 
different from that of a blue star attenuated by 
interstellar extinction law. At optical wavelengths 
the presence of scattered light - a blue component 
- may be important and at NIR wavelengths hot 
grains could have a non-negligible contribution to 
the spectrum. A model taking in account all these 
components will be explored in a forthcoming pa- 
per. 

4. Discussion on polarimetric results 

In this section, we present the estimates to the 
foreground and intrinsic polarization of Hen 3- 
1475. We also discuss a probable alignment be- 
tween the intrinsic polarization of this object and 
the interstellar magnetic field. Finally, numerical 
modeling of Hen 3-1475 polarization is done as- 
suming it is produced in a circumstellar disk. 

4.1. Intrinsic polarization of Hen 3-1475 

The observed polarization in astrophysical ob- 
jects is usually the sum of two components: (i) 
the intrinsic one, produced in the object itself; (ii) 
the foreground component, which is produced in 
the interstellar medium (ISM) between the object 
and the Earth. The high polarization of Hen 3- 
1475 suggests that at least part of it may have an 
intrinsic origin because it is very unusual such a 
high value (« 8% in V band) to be produced in 
the ISM. To obtain the intrinsic component, it is 
necessary to estimate the foreground polarization 
and subtract it from the observed one. 

The foreground polarization in each band was 
estimated by the weighted average polarization of 
all stars in the field-of-view. The results are pre- 
sented in Table 1. As the line-of-sight to Hen 3- 
1475 is near the Galactic center {I = 9.8°; 6 = 
-1-5.3°), the images contain a large number of stars 
(see last column of Table 1). To illustrate the dis- 
tribution of the polarization of the field stars we 
present histograms for their moduli and position 
angles in Ic band (Fig. l.a and Fig. l.b): no- 
tice that the maximum of each histogram coincides 
with the estimated values of the foreground polar- 
ization (Table 1). The histograms show single- 
mode distributions consistent with a sample of 



stars coming from a singular population with re- 
gard to (interstellar) polarization properties. This 
may be interpreted as an evidence that the fore- 
groimd dust shares the same properties all over 
this field. The histogram in Fig. l.a also shows 
how large the polarization of Hen 3-1475 (arrow) is 
in comparison with the field stars. This reinforces 
the suspicion that this object has a large intrinsic 
component. 

Another evidence to the intrinsic character of 
the Hen 3-1475 polarization comes from its wave- 
length dependence. We have fitted both the ob- 
served Hen 3-1475 polarization and that of the 
field stars with the Serkowski law (1973): 

P(A) = P™a.exp|-i^ln2^|, (1) 

where: Pmax is the maximum value of the polar- 
ization; Xmax is the wavelength where P,nax oc- 
curs; and if is a parameter related to the curve 
width and dependent on Xmax (Serkowski, Math- 
ewson & Ford 1975; Codina-Landaberry & Ma- 
galhaes 1976; Wilking, Lebofsky, & Rieke 1982). 
In the fits, we have used K = 1.66Amax + 0.01 
(Whittet et al. 1992). 

The results are shown in Table 2 and Figure 2: 
the polarization was normalized to Pmax to allow 
an easier comparison between the curves. In this 
graph it is clear that the wavelength dependence of 
the polarization of Hen 3-1475 and the field stars 
are different. Also, Xmax occurs at a short wave- 
length for Hen 3-1475 when compared to the be- 
havior of the field stars. Such a short value would 
be quite unusual if its origin was interstellar. 

Based on the discussion above, we are confident 
that Hen 3-1475 has an intrinsic polarization pro- 
duced in its circumstellar region. It was estimated 
subtracting the average polarization of the field 
stars from that observed for Hen 3-1475. These 
values are presented in Table 1. If we estimate the 
foreground polarization from a Serkowski law fit 
with a constant position angle, the changes in the 
intrinsic polarization are negligible. Hereafter, the 
term "Hen 3-1475 polarization" refers to its intrin- 
sic value. 
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4.2. Geometry of Hen 3-1475 and its rela- 
tion with the ISM 

The Hen 3-1475 jets are ahgned at approxi- 
mately PA 135° (Riera et al. 1995; Bobrowski et 
al. 1995). The PA of the polarization, 42.1° (in V 
band), is consistent with the observed bipolar ge- 
ometry: it is perpendicular to the jets and parallel 
to the dust torus. This behavior was already ob- 
served in young stellar objects with bipolar geom- 
etry (Hodapp 1984; Sato et al. 1985). The mech- 
anism causing this polarization will be discussed 
in the Section 4.3: in this section we would like to 
discuss a possible correlation of the geometry of 
Hen 3-1475 with the orientation of the magnetic 
field of the surrounding ISM. 

Curiously, the intrinsic polarization of Hen 3- 
1475 is nearly parallel to the polarization of the 
field stars (48.3° in V band). The latter effect 
is usually attributed to the transport of starlight 
through a dichroic medium. This dichroism is 
caused by non-spherical grains aligned by the in- 
terstellar magnetic field. Considering the stan- 
dard model for grain alignment in the ISM, the 
polarization angle is parallel to the component of 
the magnetic field in the plane of the sky (Davis 
& Greenstein 1951; but see also Lazarian 2000). 
The foreground polarization direction in the field 
of Hen 3-1475 is almost parallel to the Galactic 
plane, which is at PA 33°. This result is consistent 
with those from maps of optical polarization that 
indicate that the direction of the magnetic field in 
the Galaxy follows the Galactic plane (Mathewson 
& Ford 1970; Axon & Ellis 1976). 

The alignment between Hen 3-1475 geometry 
and field stars polarization can be interpreted as: 
(i) a chance effect; (ii) the foreground polarization 
is caused by a cloud physically near Hen 3-1475 
and in this case the IS magnetic field aligning the 
grains is also present in the Hen 3-1475 environ- 
ment; (iii) the cloud producing the polarization 
of field stars is far from Hen 3-1475 but its mag- 
netic field is aligned with the Galactic plane, so 
the alignment of Hen 3-1475 may be regarded as 
relative to the Galactic plane. 

Considering objects of similar evolutionary sta- 
tus, it is not clear if there is a correlation be- 
tween the nebula axis and the Galactic plane. Re- 
cent results indicate that the symmetry axis of 
PNe are randomly distributed (Corradi, Aznar, 



& Mampaso 1998). However, this work excluded 
PNe with small scale structures as the point sym- 
metrical knots seen in Hen 3-1475, for instance. 
On the other hand, previous works had suggested 
some correlation between PNe axis and the Galac- 
tic plane (Melnick & Harwit 1975; Phillips 1997). 
Nebulae ejected by massive stars and bilateral 
supernova remnants seem to have their axis of 
symmetry along the Galactic plane (Hutsemek- 
ers 1999; Gaensler 1998). Possible mechanisms 
that could cause such a correlation are discussed 
in Melnick & Harwit (1975) and Phillips (1997). 

Another possibility is that the axis of symme- 
try of Hen 3-1475 is determined by the IS magnetic 
field. The dust disk of Hen 3-1475 is probably a 
remnant of the AGB wind. It is quite unlikely 
that the IS magnetic field could have an impor- 
tant role in causing the asymmetry of the dense 
AGB wind^. Alternatively, the connection be- 
tween an AGB stellar wind asymmetry and ISM 
might be a fossil from the star formation phase. 
Massive stars can indeed easily retain the main- 
sequence angular momentum (Garcia-Segura et al. 
2001). We briefiy present some results in the lit- 
erature concerning the alignment of young stel- 
lar objects geometry with the IS magnetic field. 
This is still a mostly open subject in which stud- 
ies using different techniques provide discordant 
results. Dyck & Lonsdale (1979) are among the 
first authors to notice a possible alignment be- 
tween the direction of the IR polarization of pro- 
tostellar sources and the IS magnetic field. Heck- 
ert & Zeilik (1981) show that IR polarization of 
YSOs tends to be parallel or perpendicular to the 
IS magnetic field. The large scale molecular out- 
flow tends to be aligned with the IS magnetic field 
(SneU, Loren, & Plambcck 1980; Hodapp 1984; 
Cohen, Rowland & Blah, 1984; Vrba et al. 1986). 
Recent sub- and millimetric polarimetric studies of 
prcstellar objects arc not yet conclusive whether 
there is (or not) a correlation between the core 
elongation/outflow and the magnetic field of the 
containing filament (Ward-Thompson et al. 2000; 
Matthews & Wilson 2000; Glen, Walker, & Young 
1999). 

^Dgani (1998) has done some calculations showing that fil- 
amentary structures seen in old (=low density) PN can be 
explained on basis of their interaction with the magnetic 
field of the ISM. 
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4.3. Modeling the circumstellar environ- 
ment of Hen 3-1475 

To model the spectral dependence of the polar- 
ization in Hen 3-1475, we have used a three di- 
mensional Monte Carlo code for radiative transfer 
in stellar envelopes (Rodrigues 1997; Rodrigues & 
Magalhaes 2000). It includes scattering and ab- 
sorption by dnst particles, but their emission is 
neglected: therefore all the photons are originated 
in the central source. As we are modeling optical 
emission this can be considered as a reasonable; 
approximation. The grains are considered spher- 
ical and homogeneous. We have used three op- 
tical constants in the models: astronomical sili- 
cate (Draine & Lee 1984 - DL94); 0-rich silicate 
from Ossenkopf, Henning, & Mathis (1992 - 092); 
and amorphous carbon from Zubko et al. (1996 - 
Z96). The phase function considered is not the 
Henyey-Greenstein approximation, but the ana- 
lytical function obtained from the Mic theory. The 
density, p{r), was assumed to depend on the dis- 
tance as r~^. 

The diameter of the Hen 3-1475 disk torus is 
about 2" (Borkoswki et al. 1997). It was sug- 
gested that the jets are collimated near the first 
pair of knots, which is located around 3" from the 
central object (Borkowski et al. 1997). Therefore 
the aperture radius used in the polarization mea- 
surements includes the disk, but not the jets. So, 
the intrinsic polarization, as obtained in this work, 
should originate mainly by processing of the cen- 
tral source light in the optically thick disk torus. 
This configuration can indeed produce a polariza- 
tion parallel to the disk as observed in objects 
with bipolar geometry (Bastien & Menard 1988; 
Hodapp 1984; Sato et al. 1985). This can be 
explained by multiple scattering in an optically 
thick disk (Bastien & Menard 1988, 1990; Whit- 
ney & Hartmann 1993; Fischer, Henning, & Yorke 
1994). We have adopted a cylindrical geometry to 
represent the Hen 3-1475 disk. 

The models have a relatively large number of 
input parameters (Table 3). In order to minimize 
the number of degrees of freedom, we have fixed a 
given parameter if it has a reasonable good pub- 
lished estimate - see in Table 3 whether a param- 
eter was considered fixed, its value and the refer- 
ence to the estimate. If Hen 3-1475 is at 8 kpc 
(Borkowski & Harrington 2001), 1" corresponds 



to a physical size of 0.04 pc (1.7 x 10^ i?©). This 
value was used for the external radius of the disk. 
Its inner radius was estimated to be about 0.01 pc 
= 4 X 10^ i?0 according to the IRAS fluxes (Bo- 
browski et al. 1995). This is consistent with sim- 
ulations of formation of PN envelopes (Mellema 
& Frank 1995). To estimate the stellar radius we 
have considered a spectral type B3 (Knapp et al. 
1995) and a luminosity of 25 000 (Borkowski 
& Harrington 2001), which results in a value of 20 
R0. In a cylindrical geometry, there is one more 
parameter: the disk height, h, which we varied. 
We may notice that the results are invariant if all 
lengths are multiplied by a given factor and if the 
optical depth, r, is kept constant. We can also fix 
the disk inclination, adopting the inclination angle 
of the bipolar outflow obtained from Borkowski & 
Harrington (2001), which is 40°. 

Before analyzing the quantitative results for 
Hen 3-1475, we discuss how each parameter af- 
fects the polarization curve. In the optically 
thick regime, the optical depth controls basically 
the polarization degree and slightly the shape of 
the spectral polarization curve. The polarization 
curve obtained by a model with a distribution of 
grain sizes has a maximum at a wavelength related 
to the average grain size. The width of that curve 
is related to the width of the grain size distribu- 
tion. In Hen 3-1475 we do not see the polariza- 
tion maximum, but it may occur at smaller wave- 
lengths than those covered by our observations. 
The cylinder height, /i, plays a role that is not so 
obvious. For the present case, when we consider 
a fixed line of sight (= inclination), the height de- 
fines whether the central source light is blocked or 
not by the envelope: the more direct starlight, the 
smaller the polarization degree. Also, a thicker 
cylinder contains more matter, which also con- 
tributes to increase the polarization. Hence, the 
cylinder height strongly affects the polarization 
degree. Concerning the Hen 3-1475 polarization 
model, a small height allows the central object to 
be seen directly, producing a polarization that is 
not high enough. We conclude that the minimum 
cylinder height that is needed to fit the polariza- 
tion is /?,„„„ w 7. X lO'''i?0. 

We have initially tried to fit the wavelength de- 
pendence of the Hen 3-1475 polarization using a 
single grain size but the width of the obtained 
curves were too narrow, so a size distribution had 
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to be introduced. Wc have used the recipe of 
Mathis, Rumpl, & Nordsieck (1977 - MRN), the 
standard model for the interstellar medium grains. 
It is represented by a power- law, a', with ad- 
justable limits, amin and Umax- 

In order to fit the Hen 3-1475 polarization data, 
we have proceeded in the following way: we fixed 
the parameters related to the envelope geometry 
and the index of the grain-size power-law distribu- 
tion, q. Then, we let the the grain size distribution 
limits and the optical depth vary in order to find 
the best model in a weighted chi-square sense. The 
best fits are shown in Table 4. 

Let's initially suppose a grain size distribution 
that has an exponent, q, of -3.5. In spite of hav- 
ing an oxygen-rich envelope (Knapp et al. 2000), 
the optical properties of amorphous carbon pro- 
vide a better fit to Hen 3-1475 polarization data 
than the silicate ones. This is illustrated in Ta- 
ble 4 (see models 1, 5 and 7) and also in Fig. 
3. This graph also shows that all models tend 
to produce polarization curves narrower than the 
observed. An improvement could be in principle 
obtained using less steep power-law grain size dis- 
tributions. Varying the q index for an amorphous 
carbon model (see models 1, 2 and 3 and also Fig. 
4), we could improve the chi-square (best fit for q 
= -3.0), but the shape of the polarization spectral 
curve was not well described yet. An improvement 
could also be achieved by decreasing the cylinder 
height (see model 4 and Fig. 5). The grain sizes of 
the models in Table 4 do not differ significantly, in 
spite of being obtained using models with different 
values of q and h: so the dust sizes are relatively 
well-constraincxi. Considering the best fit (model 
4), the polarization of Hen 3-1475 can be fitted by 
amorphous carbon grains with size in the range 
0.063 - 0.22 /im following a -3.0 power-law. We 
notice, however, that silicate is not discarded as a 
possible dust composition. A more refined model 
should be used to better constrain the chemical 
composition (see below). Table 4 shows that the 
grain sizes are moderately affected by the optical 
constants. 

As mentioned above the measured polarization 
of Hen 3-1475 is parallel to the dust torus (see 
Table 1). Fig. 6 shows the degree and the po- 
sition angle of the polarization obtained with the 
model 4. In the position angle panel a value of 90° 
corresponds to a polarization parallel to the disk: 



this means that the BVRI observations are repro- 
duced. However, the optical depth decreases to 
the red and in some point the envelope becomes 
optically thin, so the polarization could become 
perpendicular to the disk. This would produce 
a behavior similar to the Type Ic objects from 
Trammell ct al. (1994): a rotation of 90° and a de- 
creasing polarization from small (optically thick) 
to large wavelengths (optically thin). To check 
that possibility, we have calculated the JHK polar- 
ization produced by model 4, which is also shown 
in Fig. 6^. We see that the polarization decreases 
rapidly in the IR and remains parallel to the disk. 
Therefore the rotation is not observed in model 

4. The same result is seen for amorphous car- 
bon models of Carciofi, Magalhaes, & Bjorkman 
(2002). The polarization modulus in NIR provided 
by the code must be seen with caution because 
the model does not include the dust (rc)emission, 
which may be important at those wavelengths. 

The absence of angle rotation in model 4 is 
caused by the large absorption of the amorphous 
carbon. If a less absorptive material was used (as 
silicate), the rotation could be produced. So the 
polarization angle restricts the envelope optical 
depth and also the index of refraction of the dust. 
This means that JHK polarization measurements 
would be useful to better determine the grain com- 
position in the envelope of Hen 3-1475. 

A more complete study of the dust torus of 
Hen 3-1475 using a code that includes the radia- 
tive transfer and the thermal equilibrium is being 
carried on (C. V. Rodrigues & A. C. Carciofi, in 
preparation). The flux from optical to IR wave- 
lengths, as well as the polarization, will be mod- 
eled. This will put more constraints on the physi- 
cal properties of the dust envelope of Hen 3-1475. 

5. Conclusions 

In this work we analyze BVRI polarization and 
IR photometry of the proto-planetary nebula Hen 
3-1475. The main results are summarized below: 

• the optical broad band fluxes of Hen 3-1475 
are highly polarized. A large fraction of this 



The position angle in the K band obtained using a Monte 
Carlo code has a very large statistical error (see error bars) 
because the polajrization modulus is small. 
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polarization may be produced in the circum- 
stellar environment of this object; 

• the BVRI polarization position angle is per- 
pendicular to the axis of the bipolar struc- 
ture seen in optical images. As the aperture 
used to measure the polarization is small and 
therefore does not include the jets, that im- 
plies that the central disk must be optically 
thick in all the bands observed or the dust 
must be composed of a highly absorptive ma- 
terial; 

• the wavelength dependence of the polariza- 
tion of Hen 3-1475 was modeled using a 
Monte Carlo radiative transfer code. We 
could find a good fit using an index -3.0 for 
the power-law grain-size distribution, whose 
limits are 0.063 and 0.22 fim. In spite of 
the claim for an oxygen chemistry envelope 
(Knapp et al. 2000), the optical constants 
of amorphous carbon produce a better fit 
to the optical polarization than those of sili- 
cates. However, the simple model used here 
does not allow to discard silicate dust in the 
Hen 3-1475 circumstellar envelope; 

• our CCD images provide us with polarimet- 

ric data of a large number of stars in the 
Hen 3-1475 field. So we could reliably esti- 
mate the foreground polarization in this di- 
rection. Curiously, the foreground polariza- 
tion has the same direction of the intrinsic 
polarization of Hen 3-1475. This might be 
interpreted as an interplay between the IS 
magnetic field and the Hen 3-1475 geome- 
try; 

• the IR colors of Hen 3-1475 suggests an in- 
terstellar reddening in the V band of aroimd 
3.2 mag in the line of sight. The circumstel- 
lar extinction was estimated as A(V) « 9.0 
mag. 
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Fig. 1. — Histograms of the field stars polarization 
in Ic band: a. Polarization degree; b. Position 
angle. 
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Fig. 2. — Polarization data and Serkowski fits to 
the polarization of Hen 3-1475 and field stars. 
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Fig. 4. Models to the intrinsic polarization of 
Hen 3-1475 varying the index of the power-law size 
distribution. 
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Fig. 5. Models to the intrinsic polarization of 
Hen 3-1475 for different cylinder heights. 




Fig. 6. — Polarization (a) and position angle (b) 
in the BVRIJHK bands for the best fit. The error 
bars in panel b are estimates to statistical error of 
the model. 
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Table 1 

Polarization of Hen 3-1475 and of field stars 



Object 


Band 


P 




e 


a 


Number of 






(%) 


(%) 




n 


observations 


Hen 3-1475 


B 


9.12 


0.65 


43.7 


2.0 


1 


(Observed) 


V 


8.31 


0.28 


43.1 


1.0 






Rc 


7.17 


0.29 


41.3 


1.2 


1 




Ic 


5.70 


0.20 


40.2 


1.0 


1 


Hen 3-1475 


B 


7.81 


0.65 


43.7 


2.4 




(Intrinsic) 


V 


7.03 


0.28 


42.1 


1.1 






Rc 


5.95 


0.29 


40.2 


1.4 






Ic 


4.71 


0.20 


38.4 


1.2 




Field 


B 


1.313 


0.032 


43.45 


0.7 


95 


Stars 


V 


1.298 


0.009 


48.30 


0.2 


501 




Rc 


1.243 


0.008 


46.25 


0.2 


617 




Ic 


1.041 


0.006 


47.99 


0.2 


901 



'^average of the two observations 
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Table 2 

Serkowski law fits to the observed polarization of Hen 3-1475 and to the field stars (K 

= 1.66 Xmax + 0.01 - Whittet et al. 1992) 



Object 


p 

max 




^max 


a 


K 


a 






(%) 


(%) 












Hen 3-1475 


10.8 


3.3 


0.24 


0.18 


0.40 


0.30 


0.74 


Field stars 


1.328 


0.012 


0.4746 


0.0092 


0.798 


0.015 


5.4 
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Table 3 

Input parameters for the polarization model of Hen 3-1475 



Input parameter 


Fixed 


Value 


Source radius 


Yes 


20 R0* 


Internal radius of the disk, iij 


Yes 


4 X 10^ Rq^ 


External radius of the disk 


Yes 


1.7 X 10^ Rq'^ 


Height of the disk, h 


No 


4-9 xlO^ Rq 


Inclination angle between disk axis and plane of sky 


Yes 


40° 


Optical depth in V band, Ty 


No 


Fitted 


Mininum limit of the size distribution, amin 


No 


Fitted 


Maximum limit of the size distribution, Umax 


No 


Fitted 


Exponent of the size distribution, q 


No 


-3.5 to -2.5 


Refractive index 


Yes 


Literature** 



''Estimated in this work 

''From the IRAS fluxes (Bobrowski et al. 1995) 
■^From optical image (Riera et al. 1995; Bobrowski et al. 1995) 
'^From optical spectroscopy (Borkowski & Harrington 2001) 
''See text 
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Table 4 

Parameters of the models for the spectral dependence of the polarization of Hen 3-1475 



Model 


Material 


q 


/i/2 


TV 














(Ro) 




{lim) 






1 


Amorphous carbon - Z96 


-3.5 


8 X 10^ 


9.0 


0.063 


0.21 


13.0 


2 


Amorphous carbon - Z96 


-3.0 


8 X 10^ 


9.0 


0.064 


0.21 


7.5 


3 


Amorphous carbon - Z96 


-2.5 


8 X 10^ 


9.0 


0.062 


0.21 


11.6 


4 


Amorphous carbon - Z96 


-3.0 


7 X 10^ 


12.5 


0.063 


0.22 


2.5 


5 


Silicate 0-rich - 092 


-3.5 


8 X 10^ 


9.0 


0.072 


0.26 


30.7 


6 


Silicate 0-rich - 092 


-3.0 


8 X 10^ 


9.5 


0.070 


0.26 


19.2 


7 


Astronomical silicate - DL84 


-3.5 


8 X 10^ 


9.0 


0.064 


0.29 


39.9 
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